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c3 ■ Abstract 

The Monomer Basis Representation (MBR) method developed in the first paper is applied to 
CN I water dimer in order to illustrate its application and to show its validity. The calculations are done 



by using the SAPT-5st potential surface. Monomers are treated as rigid bodies. Radial coordinate 
is separated from the angular coordinates adiabatically. MBR method is used for solving the five 



Qj ' dimensional angular problem. Then, the results of the angular calculations are fit to a Morse 

^' 
O ' function to find the potential surface for the radial motion. The results show that the method 

Q ' works efficiently and accurately. 
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I. INTRODUCTION 

Water clusters have been the subject of intense research in the hterature, since they are 
important for understanding hydrogen bonding, for interpreting the many special futures 
of the structure, dynamics and energetics of bulk solid and liquid water. Besides hydro- 
gen bonding to water molecules is of general importance in many biological, chemical and 
physical systems. 

Theoretical studies of small water clusters shows that many-body terms are very im- 
portant. For water, it is known that many-body terms may account for up to 25% of the 
interaction energy of the bulk water [l{|. It has been found that |2|: the pair interaction 
energy represents 83 — 86% of the total interaction energy of {H^O)^,. For {H-2,0)a, the same 
percentage amounts to ^ 75%, for {H20)^ only to ~ 68%. It is evident that many-body 
terms are far from negligible. It has also been found that the major part of the many-body 
cooperative effect comes from three-body terms. The four-body terms represent 2% of the 
total interaction energy of the {H20)4, and less than 4% of the total interaction energy of the 
{H20)^. Therefore, characterization of three body terms is a major step for the development 
of an accurate potential surface for bulk liquid and solid water. 

Before starting to work on three-body terms, first it is necessary to have an accurate 
description of the two body terms in the potential surface. For that purpose, water dimer 



has been studied extensively both experimentally |3l-ll6| and theoretically 
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potential surfaces has been developed |25l-l33|. 



24l | , and many 



Among the theoretical studies of water dimer. Clary and co-workers were the first ones 
to perform six dimensional calculations. Firstly, Althorpe and Clary studied water dimer 
by separating the stretching coordinate form the angular coordinates adiabatically |l7|. 
This adiabatic approximation was justified later with a more exact treatment 22|. In the 
calculations, Wigner rotation functions were used as angular basis for each monomer, and the 
molecular symmetry of water dimer was fully exploited. Later, the same authors developed 
a new method named Discrete Variable Representation Iterative Secular Equation (DVR- 
ISE) which is an extension of the diagonalization truncation method by using the variation- 



perturbation theory. The method utilizes the iterative secular equation (ISE) method 34 1. 



In the calculations, functional bases were used for some of the angles and DVR bases were 



used for the other angles. Later, Gregory and Clary studied water dimer [19| with Diffusion 



Monte Carlo (DMC) method [35 



. Leforestier and co-workers were the first ones to do 



fully coupled six dimensional calculations with basis sets. They published two papers 20[|2l|. 
The first one included some erroneous data because of a mistake in the moments of inertia 
of water monomers which were being corrected in the second one. The calculations were 
done with a coupled product basis of Wigner rotation functions by using the pseudospectral 
split Hamiltonian (PSSH) formalism in which the kinetic energy terms are evaluated in the 
coupled product basis of Wigner rotation functions and the potential energy is evaluated in 
the grid basis. Later, Chen and Light also made fully coupled six dimensional calculations 
of water dimer 22| by using the sequential diagonalization truncation method in which 
angular and the radial parts of the Hamiltonian were diagonalized successively. In these 
calculations, authors used the same decoupled product basis of Wigner rotation functions 
with the first calculations of Althorpe and Clary [l7|, and the symmetry of water dimer 
was exploited fully in calculations. Results showed that the adiabatic approximation of 
the Althorpe and Clary is successful for predicting the tunneling splittings. Finally, van 
der Avoird and co-workers developed a new potential surface called SAPT-5s 23[ by using 

]. This potential surface 



the Symmetry Adapted Perturbation Theory (SAPT) [2|, 1371 . 

was tuned for predicting the vibration-rotation-tunneling levels of water dimer which led 

n 

to the development of a new potential surface called SAPT-5st [241]. The tun ed p otential 



surface describes the experimental data with near spectroscopic accuracy 23 
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39|. In 



calculations, first an equally spaced grid was taken for the stretching coordinate and the 
angular Hamiltonian was solved at 49 different grid points by using a coupled product basis 
of Wigner rotation functions. Then, a three point contracted DVR 40|] basis were used for 
the stretching coordinate. 



Since some experimental data about water trimer is available 41r|47|, many researchers 



started to work on water trimer 
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57|. However, since the size of the water trimer problem 



is too big to handle quantum mechanically, all the theoretical studies done up to now are 
based on very reduced dimensionality models, even for rigid water molecules. The method 
developed in the first paper may make it possible to study clusters as big as water trimer 
with more realistic models as it will be seen. 

In the following section, the MBR method developed in the first paper will be used 
for calculating the vibration-rotation-tunneling (VRT) spectra of water dimer in order to 
illustrate the application of the method, and to show its validity. In calculations, adiabatic 



approximation of Althorpe and Clary, PSSH formalism of Leforestier and co-workers and 
the potential surface of Groenenboom et. al. will be used. The main difference of the 
calculations here from the previous calculations is the generation of optimized bases for 
each monomer in the cluster by using the MBR method. It will be seen that the use of the 
MBR method leads to successful results with a basis which has a much smaller size than 
any of the bases used in previous studies of water dimer. 

II. STRUCTURE OF WATER DIMER 

Water clusters are highly nonrigid, even with rigid monomers. Their potential surface 
contains more than one global minimum which are equivalent by symmetry. The potential 
barriers between these minima are not high, so that the molecule can tunnel through these 
barriers which causes tunneling splittings in the vibration-rotation-tunneling (VRT) spectra 
of the molecules. Calculation of these splittings is a good test for the accuracy of the 
available potential surfaces. Since the wave functions corresponding to the splitting levels 
has considerable amplitude in the region where the tunneling occurs, a good prediction of 
the tunneling splittings shows that the potential surface used in the calculations is accurate 
not only around the minimum (which is usually the case when there is no tunneling) but 
also in the tunneling region. 

The structure of water dimer was first determined by Dyke and co-workers via rotation 
spectra 31,14]. By examining the experimental data. Dyke realized the presence of tunneling 
splittings and made a group theoretical classification of the tunneling- rotational levels J58| by 
using Permutation Inversion (PI) group theory 59|, |60|. It was shown that, the equilibrium 
structure of water dimer has a plane of symmetry and a nonlinear hydrogen bond. The 
equilibrium structure of water dimer is depicted roughly in figure HI The — distance 
reflecting the length of the H-bond is 2.952 A, and the dissociation energy of the hydrogen 
bond, De is 3.09 kcal/mole 43 1. 

By including all feasible permutation inversion operations it is possible to generate 16 
different configurations. Due to the presence of a plane of symmetry in the equilibrium 
structure, there is twofold structural degeneracy and only eight of these structures are non- 
superimposeable. There exist three distinct tunneling motions that rearrange the H-bond 



network on time scales ranging from 1 /is to 1 ps 6ll]. The tunneling motions connect eight 



FIG. 1: Correlation diagram for the rotation-tunneling states of {H20)2 for J = 0. In the figure AS, 
I and B refers to acceptor switching, interchange tunneling and bifurcation tunneling respectively. 
Levels are labeled with the irreducible representations of the Gig PI group. 
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degenerate minima on the intermolecular potential surface (IPS). These tunnehng motions 
are acceptor switching (AS), interchange tunnehng (I) and bifurcation tunnehng (B). 

Acceptor switching has the smaUest barrier for the tunnehng. In that motion the pro- 
tons of the acceptor monomer exchange their positions. This tunnehng does not break the 
hydrogen bond. Its causes sphtting of each ro-vibrational level into two levels. 

Interchange tunneling interchanges the roles of the acceptor and donor monomers. This 
motion has the second lowest barrier and it causes further splitting of the ro-vibrational 
levels. 

Bifurcation tunneling has the highest barrier. Since acceptor switching and interchange 
tunneling, together with the inversion of the dimer, resolves all the degeneracy of the water 
dimer. The effect of the bifurcation tunneling is not to split the energy levels, but just to shift 
them. This is a result of the fact that water dimer has a plane of symmetry, which causes 
two fold structural degeneracy of each ro-vibrational level. Acceptor switching, interchange 
tunneling and the inversion of the dimer already connects eight degenerate minima to each 
other. 

The splittings for J = rotational level of water dimer are shown in figure [H Each 
energy level is labeled with the irreducible representations of the Gig PI group which is the 
molecular symmetry group of water dimer. 



III. HAMILTONIAN AND THE OUTLINE OF CALCULATION STRATEGY 

The full water dimer problem is too big to handle quantum mechanically. Therefore, 
it is necessary to model the problem. In weakly bound clusters, inter-molecular and intra- 
molecular degrees of freedom have frequencies which differ by at least one order of magnitude. 
As a result of that inter-molecular and intra-molecular degrees of freedom can be separated 
adiabatically. Thus, while treating the inter-molecular degrees of freedom, the monomers 
can be considered as rigid bodies. 

The Hamiltonian for the inter-molecular motion of a nonrigid system consisting of two 



rigid polyatomic fragments can be written as 62 1 



In the equation above, R is the distance between the centers of mass of the monomers, fi is 
the reduced mass of the dimer given by 

- - — — 



where Mj is the total mass of the monomers. 

In the Hamiltonian Ki is the kinetic energy operator of the monomer, which can be 
expressed in the body fixed frame of the monomer as 

K, = Aji + Bjl + Cjf„ (3) 

in terms of the angular momentum operators around the body fixed axis of the monomer, 
which is the molecular symmetry axis in this case. For calculations, the z axis is defined to 
be the bisector of the HOH angle. The plane of the molecule is defined to be the xz plane. 
Rotational constants A,B,C are given in terms of the moments of inertia as 

A=—,B = —,C = —, (4) 

ZJa; ■^-'y ■^-'z 

where Ix, ly, Iz are the moments of inertia of the monomers around x, ?/, z axes respectively. 
Since the monomers are considered to be rigid A, B and C are constants. 

In equation ([T]), K12 defines the kinetic energy operator corresponding to the end-over-end 
rotation of the dimer. It is given by 
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FIG. 2: Body fixed frame of water dimer. 
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In the equation above J is the total angular momentum of the system given by 



■J = jl+ J2 + L, 



(6) 



where jVs are the monomer angular momentum operators, and L is the angular momentum 
operator for the end-over-end rotation of the dimer. j is defined as j = ji + J2- 

In order to evaluate the K12 term, it is necessary to express all the angular momentum 
operators in a common reference system. This reference system is the body fixed coordinate 
system of the dimer, of course. The z axis of the body fixed frame of the dimer is defined 
along the line joining the center of mass of the monomers and the y axis of the body fixed 
frame of the dimer is defined to be along the bisector of the HOH angle of the acceptor 
monomer in the equilibrium configuration; see figure El In this reference frame, K12 is 
expressed as 

1 



K 



12 



For J = 0, this equation reduces to 

1 



(Ji + 32 + '^3iz32z + J1+J2- + J2+J1- - '^Jz3z - J+3~ - J J 



K 



12 
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(Jl + 32 + '^3lz32z + J1+J2- + Jl+J2-)- 



(7) 



(8) 



In the equations above ' denotes that the operator refers to the body fixed frame of the 
dimer but not the body fixed frame of the monomers. 

In order to solve the eigenvalue problem, first the stretching coordinate will be separated 
from the angular coordinates adiabatically. This is first done by Althorpe and Clary [iTJ, 



and led to successful results. 



First, the angular Hamiltonian is written in the form 

Hang = K + /^^^ + AK + AV (9) 

where the model Hamiltonians, /i^, for the three dimensional monomer problems are given 

by 

h^ = K, + rP (10) 

in which Ki is the monomer's kinetic energy operator in equation ([T]), and V^^ is the model 
potential energy surface for the three dimensional problem, which is the rotation of the 
monomers. 

By comparing equation (Q with the Hamiltonian, equation [H it is easily seen that AK = 
Ki2 and AV terms can be identified as 

AV = V- y^o _ yo^ (11) 

The results of the angular calculations will be used to find an effective potential surface 
for the radial coordinate so that the total Hamiltonian can be written as 

Details of monomer calculations, five dimensional angular calculations and radial calcu- 
lations will be given in sections llVi IVIII and IVIIIl respectively. 

IV. GENERATION OF A MONOMER BASIS 

A. Model Potential Surface 

For solving the monomer problem, it is first necessary to chose a model potential surface. 
The choice here will be the external field which is generated by the other monomer. 

When one wants to solve eigenstates of one monomer in the field of the other one, a natural 
question arises. Since the monomers in water dimer are not in the same environment their 
calculations will result in different eigenstates. Since the dimer includes both a donor and an 
acceptor, the results of neither calculation will generate a sufficient basis for solving the five 
dimensional angular problem. In order to handle this donor acceptor asymmetry of water 



dimer, the Energy Selected Basis (ESB) method will be used |63l. |64|. 
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In the ESB method, one generates the model potential surface as a marginal potential 
in which the potential energy at a point is the minimum value of the potential energy with 
respect to all other coordinates, i.e. 

V'{qk) = V^{qk,q7'')- (13) 

Since the coordinates of the monomer which causes the external field is varied to find 
the minimum potential, the marginal potential will sample both the donor and the acceptor 
configurations. Therefore, the monomer calculation should result in a sufficient basis for 
both the donor and the acceptor monomers. 

For solving water dimer problem with energy selected basis method, it is necessary to 
generate a three dimensional marginal potential. Labeling the monomers as i and j, the 
model potential for monomer i in the field of monomer j will be 

^.°(o) = v^(o,cr^^) (14) 

where Q = {ai,/3i,'ji) and Q = {aj,/3j,jj) refer to the Euler angles of monomers which 
describe the orientation of the body fixed frame of the monomers with respect to the body 
fixed frame of the dimer. 

B. Primitive Bases 

In order to solve the problem two primitive bases will be used, one functional basis and 
one grid basis. Functional basis being small and compact will be the primary basis. Grid 
basis will be used for the evaluation of potential energy. Primitive functional basis will be 
the symmetric top basis \jkm). In this basis, j is the total angular momentum, k is the 
projection of angular momentum onto the body fixed z axis of the monomer and m is the 
projection of angular momentum on to the body fixed z axis of the dimer. Symmetric top 
basis functions can be expressed in terms of Wigner rotation functions as 



\jkm) = ^^^-^D^:,{a,(5,^), (15) 

where D-l^f,{a, 13,'-/) is the the Wigner rotation function, expressed in terms of Euler angles 
a, /3 and 7. Its functional form is given by 

DL,{a,P,^) = e— rf^^,(cos/3)e-^->. (16) 



Exact functional form and the symmetry properties of the (i^^.(cos/3) function can be found 



elsewhere [65 1 



Following Light 66l-l68| , the representation of the Hamiltonian in the symmetric top basis 



can be called Finite Basis Representation (FBR), and the representation of Hamiltonian in 
the grid basis can be called Discrete Variable Representation (PVR). The application of 
DVR methodology to non-product bases is discussed by Corey et. al. [69|, l70| for the case of 



spherical harmonics, and later it is applied to the case of symmetric top basis |20|, |21|, |71 |. 
The deviation from Light's original formulation in the case of coupled bases is that one no 
longer seeks a unitary transformation between two bases. Therefore, two representations are 
no longer equivalent. Since the FBR basis is more compact it is the primary basis used in 
the calculations. Use of the DVR basis provides a simple way for the evaluation of potential 
energy matrix elements, since the potential energy matrix is diagonal in the DVR basis, and 
the matrix elements are given by the value of the potential energy at the corresponding grid 
point. 

C. Symmetry Adaptation of Basis Functions 

As discussed in the first paper, the monomer problem should be symmetry adapted to 
the direct product group of the pure permutation group of the monomer, and the inversion 
subgroup of the dimer. If one labels the hydrogen atoms as 1 and 2, then the pure permu- 
tation group of water monomers is the group G2 = {E^{12)}. The inversion subgroup is as 
usual e = {E, E*}. In this case, the direct product group becomes isomorphic to the molec- 
ular symmetry group of free water monomers, which is C2v{M). However, the E* operation 
here refers to the inversion of the full dimer system, and not to the inversion of free water 
monomers. Character table of the group C2v{M) is given in table [H 

In order to make the symmetry adaptation of the basis functions, it is necessary to find 
the effects of the permutation inversion operations to Euler angles and to the basis functions. 
A simple way of finding the transformation properties of Euler angles is given in appendix 
[Al By using the transformation properties given in that appendix and by defining the body 
fixed frame of the monomers as in figure [3[ effect of the permutation operations to Euler 
angles can be found. The results are given in table IIVI By using these transformations and 
the symmetry properties of the Wigner rotation functions, the effects of the permutation 
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TABLE I: Character table of the C2t,(M) permutation inversion group. This group is the direct 
product of the groups G2, whose character table is given in table [III ^-nd e, whose character table is 
given in table Hill In the table, T = x(^y means that the irreducible representation T of the group 
C2v{M) is obtained as direct product of the irreducible representation x of the group G2 and the 
irreducible representation y of the group e. 



C2v{M) = G2^e E (12) E* (12)* 


Ai = A(^G 


1 


1 1 1 


A2 = A(S)U 


1 


1 -1 -1 


Bi = B<S)U 


1 


-1 -1 1 


B2 = B^G 


1 


-1 1 -1 



TABLE II: Character table of the G2 permutation group. 



G2 


E 


(12) 


A 
B 


1 
1 


1 
-1 



TABLE III: Character table of the inversion group e. 



e 


E E* 


G 

U 


1 1 

1 -1 



inversion operations to the symmetric top basis functions can also be found, the results are 
given in table IIVI 

By using the projection operators of the group C2v{M), whose character table is given 
in table HI symmetry adaptation of the basis functions can be done. The symmetry adapted 
basis functions will be in the form of 

\jkm;p) = Np {\jkm) + {-iy\fkm)) , (17) 

11 



FIG. 3: Orientation of the body fixed frame of water monomers. Molecule is in the xz plane. 
Origin of the axes is the center of mass of the water molecule. 
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TABLE IV: Transformation properties of Euler angles under the effect of the symmetry operations 
of the C2i,(M) molecular symmetry group and the effect of the permutation inversion operations 
to symmetric top basis functions. 



E 



a 13 ^ 



\jkm) 



(12) a 13 Tr + -/ {-l)''\jkm) 
E* TT-aP -7 {-l)''\jk7n) 
(12)* vr - a /3 TT — 7 \jkm) 



with the normalization constant given by, 



N„ 



1 



(18) 



2(1 + 4o^^mo) 

In equation (ITTl) . p is either or 1, and k is either even or odd depending on the symmetry. 
Table |V] gives the values of the parameter p for each symmetry level. 

Symmetry adaptation of the grid basis can also be done. Construction of a grid basis 
for Wigner rotation functions is discussed by Leforestier 
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7l|. According to his 



prescription, the grid basis that corresponds to an FBR basis of symmetric top basis is a 
uniform grid for the angles a and 7 whose range is (0, 27r), and Gauss-Legendre quadrature 
points for cos/3 whose range is (—1,1). This choice of grid basis corresponds to a direct 
product basis of plane wave DVR bases in angles a and 7 and a Gauss-Legendre DVR basis 
in cos/3. 

The construction of symmetry adapted DVR functions is discussed by Light and Car- 
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TABLE V: Partitioning of the symmetry adapted monomer basis for the monomer calculations. 
Form of the basis functions is given in equation (J17p . 



Representation 


P 


k 


Ai 





even 


A2 


1 


even 


Bi 





odd 


B2 


1 


odd 



rington [] 



There are two approaches. First approach is to find the symmetry adapted 



DVR functions by diagonahzing a symmetric function of the coordinate in an FBR basis 
which has the desired symmetry properties [72[. Second approach is to find the DVR ba- 
sis functions by first using the usual procedure and then to obtain the symmetry adapted 
DVR functions as hnear combinations of primitive DVR functions by using the well known 
projection operators technique. This method is first used by Carrington [73|, |7^. If the 
application of symmetry operations mixes different coordinates, then the first approach is 
not applicable. It is also not applicable at all, if the primitive basis is not a direct product 
basis. 

Since the symmetric top basis is a coupled basis, the first method mentioned above is 
not applicable. Therefore, the symmetry adapted DVR functions should be constructed by 
symmetry adapting the primitive DVR functions. 

Before discussing the symmetry adaptation of the grid basis, let's define |aj), \(3j), \jk), 
as the basis functions which are localized around the points Oj, f3j and 7^ for the angles a, 
(3, and 7 respectively. Then, 

\ail3j-fk) = |ai)|/3j)|7fc), (19) 

represents a direct product basis function in the grid basis. The transformation properties of 
the Euler angles given in table [IV] requires that if at and 7^ are grid points than vr — a^ = Oj, 
TT + 7fc = 7fc, — 7fc = Tfc and vr — 7^ = %*, should also be grid points. 

Since the group C2v{M) is the direct product of the group G2 whose character table 
is given in table [III and the inversion group e whose character table is given in table lllli 
symmetry adaptation of the basis functions can be done in two steps by first symmetry 



13 



adapting to the irreducible representations of the group G2 and then symmetry adapting to 
the irreducible representations of the group e. 

Effect of the (12) operation to the Euler angles is given in table HVl With this transfor- 
mation DVR basis functions symmetry adapted to the G2 permutation group can be written 
as 

l«i/3,7fc; ^) = ^ (l«*/3i7fc) + (-l)>./3,7fc)) • (20) 

In this equation, the parameter s depends on the specific symmetry such that (—1)'^ becomes 
the character of the (12) operation in the irreducible representation that the basis function 
belongs to . The symmetry adaptation to the G2 permutation group reduces the range of 
the angle 7 by a half so that the range of angle 7 becomes 

< 7fc < TT. (21) 

Basis Functions can also be symmetry adapted to the inversion subgroup of the molecular 
symmetry group of the water dimer. In this case, symmetry adapted basis functions will be 
in the form of 

|a^/3j7fc; si) = Ni (|a,/3j7fc; s) + {-iy\a,/3j%; s)) , (22) 

where the normalization constant is given by, 

Ni = ^ ^. (23) 

V2(l + SaiCiiS-yi^^-J 

The symmetry adaptation to the inversion subgroup of the molecular symmetry group of 
the dimer reduces the range of a by a half so that the range of a becomes 

TT Sir , ,. 

2 < ". < y (24) 

The symmetry adapted DVR functions can be obtained for each irreducible representation 
of the group C2v{M) by taking different combinations of the parameters s and /. The values 
of the parameters s and / are given in table IVII for each irreducible representation. 

Before closing this section a few points about the partitioning of the grid basis should be 
mentioned. The grid points in /3 are obtained as the DVR points of Gauss-Legendre DVR 

uq . The DVR points in cos/3 are distributed 



of cos (3 by using the standard procedure [6 
symmetrically with respect to /3 = 7r/2. 
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TABLE VI: Values of the parameters s and I for the symmetry adapted grid basis. 



Representation 



s I 



A2 
Bi 
B2 





1 

1 1 

1 



In the case of the angles a and 7, grid points are evenly and periodically distributed 
between and 27r. Thus, the grid points of the angles a and 7 which are the plane wave 
DVR points are given by 

^,=.(1) (25) 

where is either a or 7, A^ is the number of DVR points and j = 0, 1, . . . , A^— 1. According to 
this prescription is always a grid point. As a result of that, the transformation properties of 
Euler angles given in table I VI 11 requires that vr should always be a grid point, too. According 
to equation ( 125|) . this is possible only if A^ is an even number. Therefore, A^ should be an 
even number in calculations. 



D. Transformation Matrix And The Evaluation of 
Potential Energy Matrix Elements 

The use of DVR provides a simple way for the evaluation of potential energy, since the 
potential matrix is diagonal in the DVR basis. However, in the case of coupled bases FBR 
basis is more compact, so that it is the primary basis used in the calculations. Consequently, 
one needs to transform the potential matrix from the DVR basis to the FBR basis, so that 
the total Hamiltonian can be evaluated as 



^{FBR) ^ ^^{FBR) _^ y{FBR) ^ j^{FBR) _^ rj.\y{DV R)rj._ 



(26) 



The way to define the transformation matrix between the symmetric top basis and the 
corresponding grid basis is discussed by Leforestier 7l| . The transformation matrix elements 
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between the symmetric top basis, \jkm), and the grid basis are given by 



where Na and A^-^ are the number of grid points for the angles a and 7, Wi3. is the weight of 
the grid point f3j of the Gauss- Legendre quadrature for cos/3. Thus, the DVR basis functions 
are defined as 



<5.,,^.,,(a,/3,7) = J2\3km)Ti^-^^ (28) 

jkm 



27r V 2 

jkm 



Ei\/^^^*.(«'/5'^)^s:.. (29) 



v^ 2? + 1 e**^"' 6*'^''''= 

jkm V ct y J- 



(30) 



The transformation matrix can be decomposed into three one dimensional transforma- 
tions: 

where one dimensional transformations are defined as 



<" = V^v^<.(cos/3,), (32) 

pimai 



e 



ik'yt 



K. = ^^- (34) 



On the other hand, the transformation matrix is not a direct product of the one di- 
mensional transformations given above, because the transformation matrix describing the 
transformation form the basis of dl^f^{cos (3) functions to the grid basis in cos (3 depends para- 
metrically to the values of k and m. For this reason, the transformation from the dl^i^{cos (3) 
functions to the grid in cos/3 should be done first. Then, the resulting intermediate basis 
becomes decoupled, and the transformation from this intermediate basis to the three dimen- 
sional grid basis becomes a two dimensional Fourier transformation. This idea is suggested 
and used by Leforestier and co-workers in their water dimer paper [20|. 
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TABLE VII: Transformation Properties of Euler angles under the effect of tlie symmetry operations 
of the GiQ permutation inversion group. 



E ai,/3i,7i a2,/32,72 

E* TT- ai,/3i,-7i vr-a2,/32,-72 

(12) ai,/3i,7r + 7i a2,/?2,72 

(34) ai,/3i,7i a2,/32,7r + 72 

(a6)(13)(24) -a2,7r- /32,7r + 72 -ai,7r - /3i,7r + 71 

V. GENERATION OF BASES FOR THE SECOND MONOMER 

Since the bases of the second monomer wih be generated from the bases of the first 
monomer. The primitive functional basis of the second monomer will also be Wigner rotation 
functions. Thus, the primitive basis of the angular problem will be the direct product of 
two symmetric top bases: \jkm) (g) \jkm). These basis functions can be written in terms of 
Wigner rotation functions as 

\j,hm,j,hm,) = ^^^^' ^8?^^^^' ^ --P^.VK^/3i.7i)^^:.,(c:^2,/32,72) (35) 

where (a,, Pi, 7^) are the Euler angles of the monomer i with respect to the body fixed frame 
of the dimer. 

In order to find the effect of the symmetry operations to the basis functions, it is nec- 
essary to find the transformation properties of Euler angles under the action of symmetry 
operations. The effect of the symmetry operations to the Euler angles are given by Althorpe 



and Clary [17|]. However, their definition of the body fixed axis for the monomers is different 
from the one used here. In their definition the water molecule is in the zy plane, while here 
the water molecule is in the xz plane. For that reason, the transformation properties of the 
Euler angles are re-derived. The difference between the transformations they have and the 
transformations given here is only in the effect of the E* operation. Transformation prop- 
erties of Euler angles is given in table IVIIi For a discussion of the derivations see appendix 

El 

By using the symmetry properties of the Wigner rotation functions, the effect of the 
symmetry operations to primitive basis functions can be found easily. The results are sum- 
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TABLE VIII: Effects of the symmetry operations to symmetric top basis functions. 



E \j1k1m1J2k2m2) 

E* {-f'+'''\jikimd2k2rn2) 

(12) {-)''' \j1k1m1j2k2m2) 

(34) {-)^^\jikimiJ2k2m2) 

(o6)(13)(24) i-y^+^^\j2k2m2Jikimi) 



marized in table IVIIII 

As discussed in the first paper, the bases of the second monomer should be generated 
from the bases of the first monomer by using the generator of the group that contains the 
permutations of the identical monomers. For water dimer the generator of the group that 
contains the permutations of the monomers is the operation (a6)(13)(24). The effect of this 
operation to the Euler angles and to the primitive bases are given in tables IVIII and IVIIII 
respectively. 

If the i^^ basis function belonging to the irreducible representation T of the group 6*2^ (M) 
of the monomer a has the form 

rt^ =J2cl\jikimi); (36) 

then, the corresponding basis function belonging to the representation T of the group C2v{M) 
of monomer b will have the form 

rf = (a6)(13)(24)rf) = J2i-^y<\jiki^i)- (37) 

I 

Although, the equations above are written for primitive basis functions, it is obvious that 
they apply to the symmetry adapted basis functions given in table |Vl 

It is also necessary to generate a grid basis for the monomer b, from the grid basis of 
monomer a. This will be done in the same way that the spectral basis of monomer b is 
generated from the spectral basis of monomer a. Thus, in order to generate a grid basis for 
the monomer b the operation (a6)(13)(24) should be applied to the grid basis functions of 
the monomer a. 
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TABLE IX: This table shows which monomer bases should be combined for obtaining bases for the 
water dimer calculations with the group Giq. In the table, labels of the irreducible representations 
are used to imply basis functions belonging to that symmetry. For an explanation of how to obtain 
mutually orthogonal basis for the doubly degenerate levels (i.e. E^ , Ey) see reference 



M 



Gi6 Bases Gig Bases 

Af {Ai (g) Ai) © {A2 X A2) A^ {Ai (g) A2) © {A2 © Ai) 
A+ {Bi © Bi) © {B2 © B2) A^ {Bi © B2) © {B2 © Bi) 
B+ {Ai © Ai) © {A2 © A2) B^ {Ai © A2) © {A2 © Ai) 
B+ {Bi © Bi) © {B2 © B2) B^ (Bi © B2) © {B2 © 5i) 
E+ {Ai © ^2) © {A2 © Si) ^- (^1 © Bi) © (A2 © B2) 
E+ {B2 © Ai) © {Bi © ^2) S- {Bi © yli) © (^2 ® ^2) 



,» /qW ^» 



If \al , (3 J ,'yl, ') is a grid basis function of monomer a, then the corresponding grid basis 
function, \a\ , /3] , 7^ ), of the monomer b will be 



\a 



I ' 



Pr^in = (a6)(13)(24)|af\/3f\7rV (38) 



The transformation properties of the Euler angles, given in table I VI II says that the 
transformation properties of Euler angles under the operation of the permutation operation 
(a6)(13)(24) is given by 

(a6)(13)(24)(ai, /3i, 71, ^2, ^2, 72) = (-"2, tt - /^z, vr + 72, -en, tt - /3i, tt + 71). (39) 

Therefore, if a grid basis function of the monomer a is localized around the point (oj, /3j, 7j), 
then the corresponding grid basis function for the monomer b will be localized around the 
point —ai, IT — /SijTT + '-fi. The effect of the permutation operation (a6)(13)(24) to the Euler 
angles of the monomers is to relabel the angles so that they belong to monomer b, and to 
change the point the basis function is localized. Thus, this operation also mixes the order 
of the DVR functions. 
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VI. COMBINING MONOMER BASES 

In the first paper, the way to combine monomer bases for water dimer is discussed and the 
bases that should be used for each symmetry are found. The results are given in table IIXI 
Thus, after choosing the right bases for each symmetry, the basis functions can be symmetry 
adapted to an irreducible representation F of the group Gie which is the molecular symmetry 
group of the water dimer by application of the projection operator 

P=l{E + x^[(a6)(13)(24)]*(a6)(13)(24)), (40) 

where x'" [('3^^) (13) (24)] is the character of the permutation operation (a6)(13)(24) in the 
irreducible representation F. These characters can be found in the character table of the 
group Gi6 which is given in table |Xl 

In order to illustrate symmetry adaptation of the basis functions, consider the irreducible 
representation Af. Table HXl says that the bases that should be used in the calculations 
should be {Ai (g) Ai) © {A2 ® A2). Therefore, the basis that will be used in the calculations 
of the Af levels is formed by taking a direct product of the Ai basis for monomer a and 
the Ai basis for monomer b and combining this with the direct product of the A2 basis for 
monomer a and the A2 basis for monomer b. If A^" is the i^^ basis function which has the 
Ai symmetry for the monomer a, and A[^ , ^3° , ^2/ are defined similarly; then, the basis 
functions of Af calculation before the symmetry adaptation to the permutation group G2 



will be in the form of either 



or 



<'' = 4?4'. («) 



iff' = ^£'4'- m 

After the symmetry adaptation to the permutation group Gg ! the form of the basis func- 
tions will be either 

^f^ = A^Afj+A^A^, (43) 

or 

^lf^ = A^4!+^^- (44) 

Thus, the symmetry adaptation reduces the sizes of both the (^i i^ A\ ) basis and the 
{A2 ® ^2 ) basis to the half of their original size. The sign between the two terms on the 
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TABLE X: Character table of the Giq PI group. This group is isomorphic to the D^fi point group. 
In the table, T = x f^ y means that the irreducible representation T of the group Giq is the direct 
product of the irreducible representation x of the group Gg, which is the pure permutation subgroup 
of the water dimer, and the irreducible representation y of the inversion group e, whose charter 



table is given in table Hill This character table is taken from the reference [58(. The correlations 
between the irreducible representations of the group Giq with the irreducible representations of its 
subgroups Gg and e are added by the author. 

(12) (a6)(13)(24) (a6)(1324) (12)* (a6)(13)(24)* (a6)(1324)* 

GiQ = G8(^e E (34) (a6)(14)(23) (a6)(1423) (12)(34) E* (34)* (a6)(14)(23)* (a6)(1423)* (12)(34)* 



At = Ai(S)G 


1 


1 


1 










1 






A^ = A2(^G 


1 


-1 


-1 








-1 


-1 






B+ = Bi®G 


1 


1 


-1 


-1 








-1 


-1 




B+ = B2®G 


1 


-1 


1 


-1 






-1 


1 


-1 




E+ = E®G 


2 











-2 












-2 


A-{ = Ai®U 


1 


1 


1 






-1 


-1 


-1 


-1 


-1 


A^ = A2®U 


1 


-1 


-1 






-1 




1 


-1 


-1 


B- = Bi®U 


1 


1 


-1 


-1 




-1 


-1 


1 




-1 


B^ =B2<S)U 


1 


-1 


1 


-1 




-1 




-1 




— 1 


E~ = E®U 


2 











-2 


-2 











2 



right hand side of the equation is plus because the character of the (a6)(13)(24) operation 
in the A^ representation is +1. 

Equations fHT]) and f H2|) apply to the B^ irreducible representation as well as it can be 
seen from table IIX[ However, symmetry adapted combinations will be different of course. 
Since the character of the (a6)(13)(24) operation for the Bf operation is —1, symmetry 
adapted basis functions will be in the form of either 

*lf' = 4?'»S'-'»S?'»K'. (45) 
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or 

<'' = 4?<-4^£'. (46) 

When the same procedure is apphed to the A^ representation, the form of the basis 
functions before the symmetry adaptation becomes either 

^5f^^ = i?!?<, (47) 



or 



,(^2") _ r(«) rW 



^'' = B^B^1- (48) 

After the symmetry adaptation one gets 

^^^ = B^B^-B^B^. (49) 

Thus, in this case symmetry adaptation will mix the two different parts. So that one can 
take either the first product basis (i?|" ® B2 ) or the second product basis (-B2 ® Bi") 
and symmetry adapt it to get the full basis. Since symmetry adaptation mixes these two 
different product bases with each other, only one of them is sufficient to make a symmetry 
adapted calculation. 

Symmetry adapted bases for all of the other symmetries can be constructed similarly. 
It should also be noted that since the character of the operation (a&)(13)(24) is for the 
doubly degenerate levels, the projection operator given in equation (HUj) is nothing but 
the identity operation. Therefore, once the basis functions are formed, they are already 
symmetry adapted as they are and there is no further symmetry adaptation. In table 
llXt two separate bases are shown for the doubly degenerate levels which are labeled with 
subscripts x and y. These bases are orthogonal to each other and they do not mix with each 
other. In order to solve the eigenvalue problem for the doubly degenerate levels either bases 
can be used. They will have the same spectrum. The fact that there are two different bases 
for a doubly degenerate level that do not mix with each other but still give the same set of 
eigenvalues is the physical explanation of the double degeneracy, of course. 

VII. ANGULAR CALCULATIONS 

In order to solve the eigenvalue problem for five dimensional angular problem. It is 
necessary to evaluate the matrix elements of the angular Hamiltonian given in equation 
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(Q. Since the contracted basis functions of the monomers are already the eigenstates of 
the model Hamiltonians, their evaluation is easy. The term AK in that equation can be 
evaluated in the primitive basis of the monomers easily, and then can be transformed to 
the contracted bases of the monomers by using the transformation matrix which is obtained 
by solving the eigenstates of the model Hamiltonian's of the monomers. The term AV, 
can be evaluated easily in the grid basis which is a tensor product of the monomer grid 
bases, and then can be transformed to the contracted bases of monomers in two steps first 
by transforming from grid basis to the primitive functional basis and then transforming 
from the primitive functional basis to the angular basis which is the tensor product of the 
contracted bases of the monomers. 

VIII. RADIAL CALCULATION 

Once the angular problem is solved at several fixed R values, the eigenvalues for the full 
problem can be found by fitting the results of the angular calculations to a Morse function 
and solving for the eigenvalues. The Hamiltonian for this one dimensional problem becomes 

^--2)111? + "'«)■ («»' 

where V{R) is the Morse function which approximates the eigenvalues of the angular calcu- 
lations at the given R values. Since the Morse potential, given in the form, 

V{r) = D(e-2°('-^o) _ 2e-"("-'^°)), (51) 

includes three parameters, D, a, tq; it is sufficient to solve the angular problem at three 
different R values. If only three points are used, there exist a unique function which fits to 
the given data. If the calculation is done at more than three points, then the Morse function 
which fits to the data can be found by making a least squares fit. 

When the Morse fit to the function is done, the eigenvalues of the Hamiltonian, given in 
equation (1501) . can be found easily; since for J = 0, the eigenvalues are analytic. According 
to Landau and Lifschitz, the eigenvalue of the n}"^ level is given by |77| 

a , 1. 



where n takes integral values from zero to the greatest value for which the expression in the 
parentheses is positive. 
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IX. DETAILS OF CALCULATIONS 



The calculations are done by using the SAPT-5st potential surface developed by Groe- 



nenbooni et. al. 
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24j ]. The source code of this potential surface was made available 



to public by Groenenboom et. al. [39| as an EPAPS document with the document num- 
ber EPAPS:E-PRLTAO-84-060018. The source code can be obtained via ftp from the site 
ftp.aip.org under the directory /epaps/ . The mass of H2O is taken as 18.010560 and the 
moments of inertia of water monomers are taken as A = 27.8806cm~^, B = 14.5216cm~^ 
and C = 9.2778cm~^. These values are the same with the values that are used in the original 
calculations of Groenenboom et. al. 

While doing the calculations primitive functional basis of the monomers are taken as 
symmetric top basis with j < 10 and m < 8. Before the symmetry adaptation this corre- 
sponds to a basis size of ~ 1650. The number of the grid points in a and 7 are set to 26 
and the number of grid points in beta are set to 15 before symmetry adaptation. All of the 
calculations are done for J = 0. 

In monomer calculations, the spectral basis is fully symmetry adapted and the grid basis 
is symmetry adapted to the permutation of the protons but not to the inversion symmetry. 
The matrix representing the Hamiltonian operator in the symmetry adapted symmetric top 
basis is stored in memory, and the diagonalization is done directly. 

While doing the angular calculations, the angular basis which is obtained as a tensor 
product of the contracted bases of the monomers is not symmetry adapted to the full sym- 



metry of the water dimer. 
Lanczos (SAL) algorithm 



nstead, the calculations are done by using the Symmetry Adapted 



78| . The use of the SAL algorithm allows one to diagonalize more 
than one symmetries at once. In the case of water dimer problem being considered here, 
SAL method made it possible to solve for the eigenvalues of the Af and B^ levels together, 
and also A^ and -B^ levels together. This results from the fact that the angular basis of 
the A^ and B^ levels and similarly the angular basis of the A2 and i?^ levels are the same 
before symmetry adaptation (see table HXl) . In the case of doubly degenerate levels calcu- 
lations should be done separately for each level since the bases of double degenerate levels 
are unique to themselves. However, use of the SAL algorithm still makes the calculations 
faster since in the SAL algorithm projection operators are used to get the symmetry adapted 
eigenf unctions. 
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TABLE XI: A comparison of the results of MBR calculations with the results of Groenenboom 
et. al. The results of Groenenboom et. al. are taken from the table III of the reference 24l |. 



The results in the table are in cm~^. MBR results are obtained by using 100 basis functions per 
monomer. In the table a is the splitting due to acceptor tunneling; i\ and i2 are the splittings 
between A^ / B^ and A2 /B2 levels due to interchange tunneling. 



Symmetry Groenenboom et. al. MBR 


At 


-1076.8643 


-1075.2116 


E+ 


-1076.4312 


-1074.8698 


Bt 


-1076.1419 


-1074.4688 


A^ 


-1065.6333 


-1063.0106 


E- 


-1065.2540 


-1062.6818 


B2 


-1064.9825 


-1062.3926 


a 


11.19 


12.14 


h 


0.722 


0.743 


i2 


0.651 


0.618 



Angular calculations are done at three different fixed R values which are 5.38a.u., 5.53a.u. 
and 5.68a.u. The ground states eigenvalues are used to define a potential surface for the 
stretching motion. The potential surface of the stretching motion is found by making a 



nonlinear fit to the Morse function by using Newton's algorithm [79 1. 

In order to converge the results it was necessary to use 100 basis functions per monomer 
for the angular calculations. 

X. RESULTS AND DISCUSSIONS 

A comparison of the MBR results with the original calculations of Groenenboom et. al. 
is given in table IXII In the table, ii is the tunneling splitting due to interchange tunneling 
between the A^ and B^ levels which is calculated as 

H = E{Bt)-E{At), (53) 



25 



FIG. 4: Convergence of the results of the MBR calculations with the number of angular basis 
functions per monomer. 
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FIG. 5: A comparison of the results of MBR calculations (lower numbers) with the original cal- 
culations of Groenenboom et. al. (middle numbers) and the experimental data (upper numbers). 
Experimental data is not available for the acceptor switching. 




where E{x) denotes the energy of the level x; Z2 is the interchange tunneling between the 



A2 and B2 levels which is calculated as 



i2 = E{B2)-E{A2)-, 



(54) 
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and a is the tunneling splitting due to acceptor switching which is calculated as 

(55) 



E{At) + E{Bt) E{A^) + E{B, 



As it can be seen from the table the results are in good agreement with each other. Es- 
pecially, the tunneling splittings are in very good agreement. From the table, it can be 
seen that the MBR calculation leads to eigenvalues which are higher than the results of 
Groenenboom et. al. This can be attributed to the fact that the stretching coordinate is 
treated in different ways in two calculations. In the calculations of Groenenboom et. al. the 



stretching coordinated is handled with a DVR grid with 49 equally spaced points [2J|. On 
the other hand, in the MBR calculations stretching coordinate is separated from the angular 
coordinates adiabatically. Therefore, because of this difference the calculations of Groenen- 
boom et. al. are less approximate than the MBR calculations. The adiabatic separation 
of the stretching coordinate from the angular coordinates was first done by Althorpe and 
Clary [17|. Their calculations are done with a different potential surface. A comparison of 
their results with other less approximate calculations which is done with the same potential 
surface is available [22|. The comparisons shows that adiabatic approximation is successful 
in predicting the tunneling splittings. The MBR results also shows that it is possible to get 
good results with adiabatic approximation. 

Convergence of the results with the number of basis functions per monomer that is used 
in the angular calculations is shown in figure HI The comparisons of the MBR results with 
the results of Groenenboom et. al. and the experimental data |lll . Il2l . Il5| is also shown 



schematically in figure [51 

From the comparison of the MBR results with the original results, it can be said that the 
MBR method gives good results. Since the number of the optimized basis functions that are 
used for each monomer (100) is much more smaller than the number of the primitive basis 
functions (~ 1650), it can also be said that the method is efficient. 

If the calculations were done with the same primitive bases but without generation of 
any optimized bases, it would be possible to decrease the size of the basis by a factor of 16 
with the help of standard symmetry adaptation procedures since the molecular symmetry 
group of water dimer is of the order of 16. On the other hand, the use of the MBR method 
decreases the size of the basis of a single monomer with almost the same factor so that the 
size of the cluster basis becomes about 16 times smaller than what one would be able to 
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achieve with standard symmetry adaptation procedures. 

The calculations given here can be improved in several ways. Firstly, the stretching 
coordinate can be treated more accurately. This can be achieved either by using more 
points to find the Morse potential or using more exact ways to handle it as Groenenboom 
et. al. have done or by using sequential diagonalization truncation schemes. However, this 
does not really seem to be necessary since the results are quite successful. Secondly, from 
figure m it can be seen that the convergence of the results is not uniform. The changes in 
the results when the number of angular basis functions per monomer is increased from 50 to 
60 are greater than the changes in the results when the number of angular basis functions 
per monomer is increased from 40 to 50. This shows that simply taking the states with the 
lowest energies as a contracted basis is not a good idea. It might be possible to devise better 
strategies while forming the contracted bases in order to obtain the best possible contracted 
basis. Although, this does not seem to be a big problem for a six dimensional system, it 
might be important when one wants to study higher dimensional systems. 

XI. CONCLUSIONS 

Application of the MBR method has been illustrated by calculating the VRT spectra of 



water dimer by using the SAPT-5st potential surface of Groenenboom et. al. [2^. The 
calculations are done by using Wigner rotation functions as primitive bases. The use of 
the MBR method made it possible to decrease size of monomer bases by a factor of ^ 16. 
The results of the calculations are in good agreement with both the original calculations of 
Groenenboom et. al. and also with the experimental results. A detailed discussion of the 
results can be found in section |Xl 

Because of its efficiency, the MBR method can be used for studies of clusters bigger than 
dimers. Thus, it can be used for studying the many-body terms and for deriving accurate 
potential surfaces. The results of calculations in this paper are especially encouraging for 
a study of water trimer with a pairwise potential surface. For example: consider the nine 
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In that model, all of the 



dimensional Hamiltonian derived by van der Avoird et. al. 
monomers are allowed to rotate around their center of masses but the center of masses of 
monomers are fixed in space. If the trimer calculation would require about the same number 
of contracted basis functions for a monomer, then a study of the nine dimensional angular 
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problem of water trimer would require 10^ basis functions. Although, a problem of that size 
can be handled with iterative methods, it is still quite big. Nevertheless, it is reasonable to 
expect that the trimer problem can be solved with less number of optimized basis functions 
per monomer. Firstly, water trimer is much more symmetric than water dimer. Secondly, 
a study of water trimer with a pairwise potential surface will have much deeper potential 
wells. Because of the importance of the three-body terms in the potential surface of water 
trimer, a calculation with a pairwise potential surface cannot give an accurate estimation 
of the experimental data. Nevertheless, this should be a big step towards the derivation of 
three-body terms in the potential surface of bulk liquid and solid water. A qualitative model 
for a possible application of the MBR method to water trimer is already given in the first 
paper. Its implementation remains as a future work. 

Appendix A: A Simple Way to Find The Transformation Properties of Euler Angles 

The transformation properties of the Euler angles of a body fixed frame under rotations 
is well know if the Euler angles are defined with respect to a space fixed frame, which does 
not move under the effect of any symmetry operation. These transformations are given 



in the book Molecular Symmetry and Spectroscopy, written by Bunker and Jensen [59|. 
For convenience of the reader, they are given in table IXIII However, the case of the small 
clusters is different, because the Euler angles of the body fixed frames of the monomers 
are defined with respect to the body fixed frame of the cluster, and not with respect to a 
space fixed frame. Since the body fixed frame of the cluster is not fixed in space, it also 
rotates with the effects of permutation inversion operations. Therefore, the application of 
the transformations given in table IXiT] is sufficient only if the body fixed frame of the cluster 
does not move. However, they will give wrong answers if the body fixed frame of the cluster 
moves, too. For example, in the case of water dimer, transformations given in table IXIII 
will give the right answer for the operation of (12). However, they will not work for the 
operations E* and (a6)(13)(24), because these operations rotate the body fixed frame of the 
dimer around its x axis by vr radians. 

To the best of author's knowledge there is not any easy way of finding the transformations 
in these cases available in the literature. In the following lines, an easy way of finding these 
transformations will be discussed. 
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591 ] . Coordinates 



TABLE XII: Transformation properties of Euler Angles. This table is taken from page 266 of 
the book Molecular Symmetry and Spectroscopy written by Bunker and Jensen 
are relabeled to make it consistent with notation used here. In the table, i?T is a rotation of the 
molecule fixed (x, y, z) axes through vr radians about an axis in the xy plane making an angle 4> 
with the X axis {(/) is measured in the right hand sense about the x axis), and R^ is a rotation of 
molecule fixed {x,y,z) axes through 9 radians about z axis {6 is measured in the right hand sense 
about the z axis). 



Hi n^ 



(3 TT- f3 13 

a a + TT a 

7 27r — 20 — 7 7 + I 



In order to find the changes in the Euler angles, the process will be divided into two steps. 
In the first step the orientation of the body fixed frame of the monomers will be kept fixed 
with respect to a space fixed frame. However, the body fixed frame of the cluster will be 
rotated together with the body, of course. In the second step, the rotation of the monomer 
frames will be done. 

The transformations in the second step will be the same with the ones given in table 
Kill because the cluster frame does not move in the second step. Therefore, if one finds the 
transformations for the first step, then the overall transformation can be found by applying 
two transformations sequentially. 

The transformation properties of the first step can be found easily by realizing that the 
only difference in the first step is that the roles of the two frames are interchanged. This 
time, it is the monomer frame which acts like the space fixed frame since it doesn't move; 
and it is the cluster frame that rotates. If (a,/3,7) are the Euler angles of the monomer 
frame defined with respect to the cluster frame, one can equally say that (—7, — /3, —a) are 
the Euler angles of the cluster frame defined with respect to the monomer frame. Thus, 
the Euler angles of the cluster frame with respect to the monomer frame after the rotation 
can be found by applying the transformations in table IXIIi Then, the Euler angles of the 
monomer frame with respect to the cluster frame can be found by using the same trick 
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TABLE XIII: Transformation properties of Euler Angles. Definitions of the rotation operations 
are the same with table IXlIl However, this time it is the cluster frame rotating, not the monomer 
frame. 



/3 7r-/3 /3 



a Z0 — a a 



7 TT + 7 7 



again. If {a', (3', 7') are the Euler angles of the cluster frame obtained by applying the 
transformations in table IXlTl to the Euler angles (—7, — /3, —a), then the Euler angles of the 
monomer frame will be (—7', — /?', —a'). 

Thus, by using the transformation rules given in table IXI If one can find the transformation 
properties of Euler angles of the monomer when it is only the frame of the cluster that rotates. 
The general formulas for such transformations are derived and the results are summarized 
in table IXlTTl 

Since the order of the two successive operations that are used to find the transformation 
properties of the Euler angles of the monomer is not important for the final orientation of 
the frames, the order in which one uses the transformations is not important. 

To illustrate the applications of the transformations given in tables IXIIt and IXIIH to 
molecular clusters, consider water dimer. Body fixed frames of the dimer and the monomers 
are given in figures [2] and |3l respectively. Firstly, consider the effect of E* operation to water 
dimer. This operation rotates the body fixed frame of the monomer around its y axis by vr 
radians (see figure [3]) and it also rotates the body fixed frame of the dimer around its x axis 
by IT radians (see figure [2]). By using the transformations in table IXlTl (with (f) = 7r/2), the 
effect of the monomer rotation to the Euler angles is found to be 

(a,/3,7) ^ (7r + a,7r-/3,7r-7). (Al) 

Then, by using the transformations in table IXIIII (with = 0) the effect of rotation of the 
dimer frame to the Euler angles of the monomer frame is found to be 

(tt + Q;,7r - /^,7r - 7) -)■ (tt - a,/?, -7, ) (A2) 
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where —a — tt is replaced with vr — a. This is possible because adding 2tx to the angles 
a and 7 doesn't change anything since they have a period of l-n. Therefore, the overall 
transformation becomes, 

(a,/3,7)^(7r-a,;g,-7). (A3) 

This is what is reported in table IVIII The same transformation applies to both of the 
monomers. 

In the case of the operation (a6)(13)(24), the effect of the operation to the monomer 
frames is just to relabel them, so the first transformation becomes 

(ai, /3i, 7i, as, ^2-, 72) -^ ("2, ^2-, 72, "i, A, 7i)- (A4) 

This operation also rotates the dimer frame around its x axis by tt radians, so with the 
transformations given in table IXIIIt one gets 

("2, /32, 72, ttl, /^l, 7l) -^ (-"2, TT - /^a, Vr + 72, -«!, TT - /3i, TT + 7i). (A5) 

Therefore, the overall effect of the (a6)(13)(24) operation to the Euler angles of the monomers 
become 

(«!, /3i, 7i, ^2, /32, 72) -^ (-«2, vr - /32, VT + 72, -ai, vr - /3i, vr + 71). (A6) 

The symmetry operations (12) and (34) do not have any effect on the body fixed frame 
of the dimer. Therefore, their effect to the Euler angles of the monomers can be found by 
using the transformations given in table IXIIi The results will be the ones that are reported 
in table IVm 

Appendix B: Real Symmetric Top Basis 

The fact that the symmetric top basis is complex, makes the calculations more demanding, 
since complex numbers occupies twice the memory real numbers occupies. For this reason, 
it is advantageous to transform the symmetric top basis to a real basis that we will be called 
real symmetric top basis. The way to generate a real basis from a complex basis is obvious: 
sum and differentiate with its complex conjugate. By using the symmetry properties of the 



Wigner rotation functions 65|, it can be shown easily that 

^^\(«,/9,7) = (-l)"^^'^i^-fc(«,/9,7). (Bl) 
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From the equation above, it follows that Ijkm)* = (— l)"*~'^|jfcm). This leads to the 
result that the linear combinations \jkm) ± \jkrh) are either pure real or pure imaginary. 
In the case of pure imaginary functions, the imaginary number i can be omitted since it is 
just a phase factor. After doing some algebra, following functions are obtained as the real 
symmetric top basis. 

For m — k is odd: 

iJkm) + \jkfh) -\/2j + ] 



Za/2 + 25ko5mO 2-\/l + 5kQ5mQ 



sin{ma + k'j)dl^i^{cos/3), (B2) 



\jkm) - \jkm) ^/2j + 1 • 

7n = o /I I A A cos{ma + k-f)dl^{cos/3). (B3) 

For ?T2 — fc is even: 



^7^^ = ;r7f=^^^sm(ma + /c7)ci^^;^(cos/3), (B4) 

«V2 2Vl + ^fcO^mO 

— , ^^ = — =cos(ma + fc7)a:^^(cosp). (B5) 

V2 + 25fco^mO 2V 1 + ^fcO^mO 

When the real symmetric top basis functions are used as a basis, the transformation 
matrix elements from this basis to the grid basis also becomes real. Thus, by using equation 
( l27l) . transformation matrix elements for the basis functions given by equations ( 1B2[) and 
(IB4p becomes 

and the transformation matrix elements for the basis functions given in equations (1B3P and 
( JBSP becomes 

T^jfem v/2j + 1 cos(ma^ + k'jk) j 

Appendix C: Kinetic Energy Matrix Elements in Water Dimer Calculations 

In this appendix, the matrix elements are given for the complex symmetric top basis 
functions. Application of the formulas given here to the real symmetric top basis functions 
is straightforward. 

If \jikimiJ2k2'm2) = \jikimi)\J2k2'm2); then, the nonzero matrix elements of Ki is given 
by 

{jikimiJ2k2m2\Ki\jikimiJ2k2m2) = i — - — j ji{ji + 1) + ( C — j kf, (CI) 

33 



and 

75 (^ 

{jikimiJ2k2m2\Ki\jiki ± 2miJ2k2m2) 



4 
xVii(ii + l)-(fci±2)(A;i±l) 

X y^j^iJ^TTf^^lhiWi- (C2) 

Similar expressions can be derived for K2. For J = 0, nonzero matrix elements of Ku 
are given by 

{jikimiJ2k2m2\Ki2\jikimiJ2k2m2) = ■^— ^[ji(ji + 1) + J2(j2 + 1) + 2mim2]6m.^^-m.^, (C3) 
and 

xC±^/^„_^„ (C4) 



where 



c,t = Vj(j + i)±^(^ + i)- (C5) 
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